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The widely accepted dark matter hypothesis offers a seductive solution to missing mass prob- 
lems (galaxies, clusters of galaxies, gravitational collapse in structure formation,...). However 
the physical nature of the Dark Matter itself is still unknown. Alternatively, it has been pro- 
, posed that apparent dynamical evidence of dark matter is due to a modification of Newton's 

law of gravitation. Here we revisit the Modified Newtonian Dynamics (MOND) theories at the 
scale of galaxy clusters. Using hydrodynamical simulations, we derived quantities such as the 
density and the temperature of the ICM. We compared those MOND simulated predictions 
5^ ' to high quality X-ray density and temperature profiles observed down to ~0.5 the virial ra- 

dius. If the density profiles seems in acceptable agreement, the simulated temperature show a 
constant increase with the radius whereas the observed profiles show a flat to a mild decrease 
shape down to ~0.5 Rvirial- We also computed the dynamical MOND mass for 8 X-ray clus- 
ters observed with XMM-Newton. If the MOND mass helps to lower the discrepancy with 
the baryonic mass by ~20%, still ~80% of the mass in clusters is unaccounted by baryons. In 
order to solve this problem and to reconcile MOND with clusters observations, we investigated 
the possibility of an added dark baryonic component. We assumed a component of massive 
neutrinos to fill the remaining discrepancy between the observed MOND dynamical mass and 
the baryonic mass. This led us to derive a tied observational lower limit for the neutrino mass, 
m v = 1.06 eV. 



1 Introduction 

The missing mass problem in clusters of galaxies arises from the comparison of the observed 
baryonic mass with the observed dynamic mass. The baryonic mass is mainly due to the hot 
intracluster gas that is well observed in X-rays via its free-free emission. The current status of 
the observed gas fraction in clusters gives a fairly well constrained value of about 12% [see 10, 
2, for instance]. Taking into account the stellar mass, this makes the discrepancy between the 
observed dynamic mass and the observed baryonic mass larger than a factor of 7. 
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Figure 1: Scaled mass profiles of ten nearby clusters. The mass is scaled to M v , and the radius to R v , both 
values being derived from the best fitting NFW model. The solid black line corresponds to the mean scaled NFW 
profile and the two dashed lines are the associated standard deviation. Figure from [16]. 

We investigate in the following two hypothesis to solve this problem, the Dark Matter 
paradigm and the Modified Newtonian Dynamics paradigm. 

2 The Dark side of the matter 

The dark matter (DM hereafter) hypothesis appears to provide a seductive explanation of the 
mass missing problem. A new component of non-baryonic matter, insensitive to all interactions 
but gravitation, is introduced to fill the gap between the baryonic matter and the binding mass. 

[16] have measured the total mass profile over a sample of ten clusters from 0.01 H v up to 
0.5 R200 using XMM-Newton observations. Their sample has an excellent temperature coverage 
and covers an order of magnitude in mass from M v = 1.2 x 10 14 M to 1.2 x 10 15 M Q . 

They have found that the NFW profile is a good representation of the ten observed mass 
profiles, and that in all cases the isothermal sphere model (i.e a profile with a core) is rejected at 
high confidence. In other words, they confirm the cusped nature of the Dark Matter profile, as 
predicted by CDM simulations of hierarchical structure formation, over the temperature/mass 
range of the sample. The mass profile shape is close to universal, again as predicted, with a 
dispersion of less than 15% at 0.1 R„ in the scaled mass profiles. The shape is quantitatively 
consistent with theoretical predictions. The variation of the observed concentration parameters 
with mass is in line with the predictions, taking into account the measurement errors and the 
expected intrinsic scatter. Taken together, these results provide strong evidence in favour of the 
Cold Dark Matter cosmological scenario, and show that the physics of the Dark Matter collapse 
is well understood. The ten scaled profiles are presented on Fig. 1. 

However, while cosmological evidence is accumulating in favour of this scenario (see for 
instance [7, 22, 23]), it is disconcerting that the nature of the non-baryonic dark matter is 
completely unknown. Of course there are many candidates of varying degrees of detectability 
and plausibility (e.g review by [4]). 

3 Living in a MOND world 

As an alternative to dark matter scenarios, [13] proposed a modification of the Newtonian 
dynamics effective at galactic and extra-galactic scales. This modified Newtonian dynamics 



(MOND hereafter) has been notably successful in explaining the discrepancy between rotation 
and luminosity curves in spiral galaxies [11, 12], and claims other phenomenological successes 
(for a full review on MOND see [19]). The discrepancy between the baryonic mass and the 
dynamical mass in clusters of galaxies is perhaps foremost among the issues that MOND has 
yet to convincingly address. 

The first confrontation of X-ray observations of clusters with MOND [9] emphasised the 
difficulties faced by MOND in passing the cluster test. The problem was revisited by [17, 
18] and ended in a remaining discrepancy of a factor of 2-3 between the baryonic and the 
MOND masses. More recently, [1] discussed observational evidence for three clusters for which 
the observed discrepancy is about 1-5 within 1 Mpc and is boosted to a factor ~10 within 
the central 200 kpc, further weakening the reliability of the MOND paradigm. However, [20] 
responded with an update of his earlier work, introducing an added ad hoc dark component at 
the cluster centre to reduce the discrepancy to only a factor of 1-3 overall in the cluster. 

Some other tests have also been carried out using gravitational lensing data. They have also 
pointed out the difficulties faced by MOND at the cluster scale [8, 5]. 

4 Revisiting MOND at galaxy clusters scale 

Recently, [15] and [14] have revisited the MOND paradigm at the cluster scale. In the following 
we present a synthesis of their results. 

4-1 Structure formation in MOND 

[14] use a one dimensional hydro dynamical code to study the evolution of spherically symmetric 
perturbations in the framework of Modified Newtonian Dynamics (MOND). The code used 
evolves spherical gaseous shells in an expanding Universe by employing a MOND-type relation- 
ship between the fluctuations in the density field and the gravitational force, g. They focus on 
the evolution of initial density perturbations of the form Si ~ r~ s for < s < 3. A shell is 
initially cold and remains so until it encounters the shock formed by the earlier collapse of shells 
nearer to the centre. During the early epochs g is sufficiently large and shells move according to 
Newtonian gravity. As the physical size of the perturbation increases with time, g gets smaller 
and the evolution eventually becomes MOND dominated. However, the density in the inner 
collapsed regions is large enough that they re-enter the Newtonian regime. The evolved gas 
temperature and density profiles tend to a universal form that is independent of the the slope, 
s, and of the initial amplitude. 

The results from the simulations, the temperature and density profiles of MOND, were 
confronted with recent X-ray observations of nearby galaxy clusters by XMM-Newton [16] and 
Chandra [24]. Among the different MOND runs, the minimum value number density at 500 
kpc in MOND is 0.07 cm~ 3 . This is substantially higher that the maximum observed number 
density of 0.0026 cm~ 3 . More importantly, the temperature profile, ~ r 0,5 in MOND, is in clear 
contrast to the observed profiles, which either flatten or show a mild decline in the outer regions. 

4-2 Observed clusters mass profiles versus MOND predictions 

From the mass profiles derived by [16] in a LCDM cosmology, [15] computed the observed ratio 
of the dynamical MOND mass, M, m to the baryonic mass : 
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Figure 2: MOND simulated temperature profiles (thick black line) compared to the individual profiles of clusters 
observed with XMM- Newton (red dashed lines) and Chandra (blue dotted lines) . Bottom panel: MOND simulated 
density profiles (thick black line) shown with respect to the observed density profiles of nearby clusters by the 
XMM-Newton satellite (red dashed lines). In both panels the profiles have been normalised according to their 
respective values at r = 500 kpc as marked radius by dotted vertical line. Figure from [14]. 

where f gas is the gas fraction, f+ is the stellar fraction (i.e the ratio of the stellar mass to 
the gas mass of the cluster), ao ~ 10~ 8 cm/s -2 is a fundamental accelaration constant in the 
MOND theory (its value is derived from the analysis of the rotation curves of galaxies in a 
MOND framework), G is the gravitational constant, is the hydrostatic mass of the cluster 
as obtained from the hydrostatic equilibrium equation. 

In a pure MOND universe this ratio should equate 1. Table 1 reports the average baryon 
fraction values (i.e fb = / S as(l + /*)) and the average M m /Mb ratios for the measured clusters at 
different radius. The computation were done at the following radii: the radii corresponding to 
the density contrasts, with respect to the critical density of the Universe at the cluster redshift 
of 5 = 1000 and 5 = 15000 (i.e 0.47 ± 0.02 R 200 and 0.10 ± 0.01 R 20 o average over the eight out 
of ten clusters from [16]). Those two radii mark the boundaries of the radial range over which 
the observational constraints are especially well tied down. 

At 0.5 R„, the ratio of the MOND mass to the baryonic mass is 4.94 ± 0.50 (10.6 ± 3.77 at 

Table 1: Ratios of the dynamical mass to the MOND mass and of the MOND mass to the baryonic mass over 8 

nearby clusters. 



Radius 


fb 


M m /M b 


<5 = 15000 


0.09 ± 0.03 


10.6 ±3.77 


5 = 1000 


0.13 ±0.02 


4.94 ±0.50 



~ 0.1 R^). This is more than a factor of two above the value derived by [18]. The evidence 
is confirmed if we only consider the hot systems. Indeed, for clusters with kT > 3.5 keV, 
M d /M m = 1.43 ± 0.08 at ~ 0.5 R 20 o- This makes the ratio of the MOND mass to the baryonic 
mass ~ 5.10 ± 0.56. Thus in all cases, within a 3d (i.e. 99% confidence) the ratio M m /Mb in a 
MOND cosmology will be greater than 3.4, making MOND unable to fully overcome the missing 
mass problem in clusters. With respect to the used sample, in a MOND framework, still about 
80% of the total mass of galaxy clusters is missing at ~ 0.5 R200 • I n other terms, this means 
that MOND just reduces the missing mass problem in clusters by about 20% (at half the virial 
radius), but does not solve it. 



4-3 Filling clusters with neutrinos 

A last alternative to rescue MOND is to invoke a non-luminous component at the centre of 
clusters, as suggested recently by [20]. This author proposed massive neutrinos, aggregating at 
the cluster scale, as candidates for this dark component. More recently, [21] also called upon 
this neutrino hypothesis, studying formation of structures in the relativistic MOND framework 
(i.e. the Bekenstein theory [3]). 

Further assuming a constant density sphere for this dark component and taking into account 
the phase space density limit for neutrinos, [20] derived an upper limit for the neutrino density af- 
ter their collapse and accretion within structure of: p v < (4.8xl0~ 24 )(m i ,/2eV) 4 (Tk ;e y) 3//2 kg/m~ 3 
We make use of this limit on the neutrino density, and we use the spectroscopic temperatures 
measured for each cluster of the previous sample between 0.1 and 0.5 R200 ( see [16])- From the 
eight clusters, it is possible to compute the needed neutrino mass to equate the missing mass 
at a given radius (i.e. M m (r) — Mb(r)) with the contribution of the massive neutrino to the 
cluster total mass. [20] hypothesis of the neutrino accretion mainly concerned the central parts 
of clusters. In our study, to explain the ~80% of missing mass in MOND, we add to extend 
the radius of the neutrino sphere down to R1000 • The minimum neutrino mass then required is 
m u > 1.74 ± 0.34 eV . This is a strongly constraining value for the neutrino mass, which makes 
the lower bound for the neutrino mass becomes ~ 1.06 eV, within a 2a limit (i.e. 95% confi- 
dence). This is barely compatible with the cosmological constraints from combined CMB+LSS 
data [22, 6]. 



5 Conclusion 

We have presented here the case of two hypothesis to solve the missing mass problem at the 
galaxy clusters scale: The Dark Matter and the MOND paradigms. 

The Dark Matter hypothesis remains a very secure and stable grounds as it solve by definition 
the missing mass problem. However, the question of the nature of the Dark matter is still one 
of the most dazzling problem in modern astrophysics. 

The MOND hypothesis, as a stand alone solution, has proven to be unable to solve the 
missing mass problem at the clusters scale. An added hot dark component is needed to rescue 
MOND, massive neutrinos for instance. As the amount of hot DM needed is up to ~80% of 
the clusters mass, this turns the MONDian cosmological framework more into a mixed DM 
cosmology. 

The results presented to this conference are extensively detailed the three following papers: 
Pointecouteau, Arnaud & Pratt (2005) [16], Pointecouteau & Silk (2005) [15] and Nusser & 
Pointecouteau (2006) [14]. 
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